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Using  Gaussian  plume  model.  Considering  affection  of  terrain,  pollutant 
decay  and  p>  ;cipilali«n  scavenging,  tlie  formula  for  esTTniatiug  the  long-term 
average  concentration  from  the  power  plant's  high  stack  was  proposed  in  this 
paper.  Hecause  the  wind,  stability  and  mixing  depth  are  calculate  d  from 
routing  meteorological  data,  we  can  estimate,  the  average  concentration  without 
any  field  observation.  It  is  a  significative  method  on  economy  and  practical 
use  for  electing  plant  location  and  assessing  air  quality. 


A  METHOD  OF  ESTIMATING  THE  LONG-TERM  AVERAGE  CONCENTRATION  FROM  THE 
POWER  PLANT'S  HIGH  STACK 


Ma  Fujian,  Wang  Yanchang,  Pan  Yunxian  and  Zhou  Chaofu. 


Using  the  Gaussian  plume  model,  this  paper 
considered  the  effects  of  terrain,  pollutant 
decay  and  precipitation  scavenging,  and  proposed 
the  formula  for  estimating  long-term  average 
ground  level  concentration  for  the  power  plant's 
higVTstack.  Since  parameters  such  as  wind, 
stability,  mixing  depth,  etc.  were  calculated 
from  routine  meteorological  data,  the  average 
concentration  could  be  estimated  without  any 
field  observation.  This  is  a  method  with  sig¬ 
nificance  in  practical  use  and  economy  for  selec¬ 
ting  a  preconstruction  plant  site  and  assessing 
air  quality .  -  - 


The  assessment  of  air  quality  and  selection  of  a  plant  site  for  a 
f oss i 1 -fueled  power  plant  both  require  estimating  the  long-term  average 
concentration.  Our  country  has  set  coal  as  the  primary  fuel  for  fossil 
fueled  power  plants.  Most  of  them  are  located  near  coal  mine  entrances 
or  in  hilly  regions,  and  usually  lack  pollution  meteorological  observa- 
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tion  data.  Therefore,  how  to  use  routine  meteorological  data  from 
nearby  weather  stations  to  estimate  the  long-term  average  concentration 
distribution  caused  by  emissions  from  fossil-fueled  power  plants  is  an 
urgent  problem  waiting  for  solutions." 

Currently,  the  popular  Gaussian  plume  model  discribes  the  rules  of 
point  source  diffusion  quite  successfully.  It  has  a  simple  form  and 
clear  physical  concept,  and  the  model  estimated  values  correspond 
fairly  well  with  the  observed  values;  therefore,  it  enjoys  extensive 
application.  When  this  model  is  used  to  estimate  pollutant  long-term 
average  concent raiton  distribution,  meteorological  parameters  such  as 
wind,  stability  and  mixing  depth  are  required  in  addition  to  data  in 
plume  rise  and • d i f f us  ion  coefficient.  Effects  of  terrain  and  precipi¬ 
tation  scavenging  should  also  be  considered  for  the  model.  Through 
numerous  experiments  and  theoretical  research,  many  formulae  and 
figures  have  been  proposed  which  can  be  used  for  selecting  proper 
atmospheric  diffusion  coefficients  and  plume  rise  [1],  Thus,  how  to 
use  routine  ground  level  meteorological  observation  data  to  estimate 
atmospheric  stability  and  mixing  depth  has  become  the  key  to  the 
application  of  the  model. 

This  paper  proposed  a  model  and  method  that,  with  given  diffusion 
coefficient  and  plume  rise,  uses  the  routine  ground  level  meteorologi¬ 
cal  data  to  estimate  the  ground  level  long-term  average  concentration 
caused  by  emissions  from  a  f oss i 1 - f ueled  power  plant's  high  stack, 
and  an  example  was  given.  Said  method  is  simple  to  use  and  is  especial 
ly  suitable  for  regions  with  only  ground  weather  stations.  It  provides 
an  effective  method  for  the  industry  to  conduct  preconstruction  site 
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selection  and  air  quality  assessment. 


I.  Establishing  The  Model 


The  well-known  Gaussian  plume  model  is: 

« p[-  4=  W'-S?  1 

where  <5  and  <S  are  horizontal  and  vertical  diffusion  coefficients 

y  x 

respectively.  H  is  the  effective  stack  height. 


In  the  atmospheric  boundary  layer,  due  to  the  thermal  and  dynamic 
action  of  its  under  side,  a  layer  with  violent  turbulence  movement  is 
formed  above  the  ground  and  is  called  the  atmospheric  mixing  layer. 

The  portion  above  the  mixing  layer  is  usually  of  stable  structure  and 
will  suppress  the  pollutant's  upward  diffusion.  The  high  stack  of  a 
fossile-f ueled  power  plant  can  reach  200  meters  making  the  effective 
stack  height  reaching  hundreds  of  meters.  Therefore,  the  model  should 
take  the  suppression  and  reflection  effects  of  the  stable  structure 
above  the  mixing  layer  into  consideration.  The  following  basic  form 


of  the  Gaussian  point  source  model  is  adopted: 

c  (.x.  y  O  :H  > 


where  L  is  the  mixing  depth,  n  is  the  number  of  reflection  -  (generally 
select  between  3  to  5). 


The  effects  of  different  meteorolog ica  1  conditions  or.  the  model 
parameters  must  be  considered  when  deriving  the  formula  for  long-term 
average  concentration.  The  wind  direction  is  divided  into  fifteen 
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In  the  above  formula,  if  m=4  which  indicates  stable  atmospheric 
condition,  then  in  Equation  (3)  n=0. 

The  formula  for  long-term  average  concentration  takes  the  effects 
of  the  mixing  layer  on  high  stack  diffusion  into  consideration.  It 
also  considers  effects  of  the  discontinuous  modification  on  the  border¬ 
ing  lines  between  the  16  wind  directions,  terrain  adjustment  proposed 
by  Egan  for  a  situation  where  the  plume's  effective  stack  height  is 
greater  than  the  terrain  height  (h^,)  ,  pollutant  decay  and  precipita¬ 
tion  scavenging. 

For  a  rough  under  side,  the  effects  of  terrain  on  the  plume's 

diffusion  can  simply  be  included  in  the  effects  of  terrain  on  the 

r  2  1 

plume's  effec-tive  stack  height.  Egan  (  1975)  proposed  to  replace  H 

with  TeH,  where  Tc  is  the  terrain  modification  factor.  When  H^hT, 

hT 

Tc=l-hT/2H,  thus  TeH-  (  H — y)  . 

L  is  the  half-life  of  pollution  decay  (in  hours).  x/u  (x  is  the 
downwind  distance  from  the  point  source,  u  is  average  wind  speed)  is 
the  time  required  (in  seconds)  for  the  pollutant  to  travel  from  the 
source  to  the  receiving  point.  These  express  the  effects  of  decay  on 
concentration  distribution  during  the  pollutant's  transport  process 
[3].  The  process  of  precipitation  scavenging  is  more  complicated,  and 
its  effects  on  long-term  average  concentration  can  be  roughly  estimated 
dv  adopting  the  average  effects  method. 

If  the  effects  of  discontinuous  modification,  terrain  modifica¬ 
tion,  pollutant  decay  and  precipitation  scavenging  are  not  considered, 
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formula  (3)  can  be  transformed  to  the  simplest  form  of  long-term 
average  concentration: 


and 


c  (x.o  .u  n >=  2  i:  i;  *tmlix,o.  a  h) 

>  *  1  m**l  /-] 


U\ml(x,  O.  J6 

2  * 


2  Q 


2  nxu,*.' 


(2 


2*7, 


-J 


(4) 


II.  Determination  Of  Parameters 


Using  formula  (3)  to  estimate  long-term  average  concentration 


requires  values  for  the  following  parameters: 

4 


xm 


H  .  ,  U.  ,  L  , 
ml  kml  m 


etc.  I_f  plume  rise  and  diffusion  coefficient  are  selected 
from  existing  research  results,  then  other  parameters  can  be  estimated 
from  routine  ground  meteorological  observation  data. 


The  most  commonly  used  method  to  divide  atmospheric  stability  is 
the  P  -  T  classf ication  method,  yet  the  cloud  height  data  required  by 
said  method  are  not  measured  by  the  ground  weather  stations  of  our 
country.  Based  on  the  characteristics  of  total  and  low  cloud  volumes 
measured  by  a  ground  weather  station  and  using  the  P  -  T  method  as  a 
base,  the  cloud  indexes  classifying  radiation  classes  are  substituted 
by  total  and  low  cloud  volumes,  and  this  is  called  the  P  -  C  method 
[4].  The  classification  results  using  this  method  correspond  consist¬ 
ently  with  those  using  the  P  -  T  meTA  od .  The  P  -  C  method  classifica¬ 
tions  are  listed  in  Table  1  and  2. 
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7)  <  4 


>  9  /  4>  i) 


Key:  (1)  Cloud  volume;  (2)  night;  (3)  Altitutde  angle  of  the  sun; 

(4)  Total  cloud  volume/Low  cloud  volume. 


Atmospheric  Stability  Classes 


)lO*fc,R 

&(*/#> 


2  — 2 , 9_ 


(2)  * 


&  »  «  *n 


A  ;  A-B 


3  -  4. 


B  1  B-C 


C-D 


Key:  (!)  wind  Speed  at  a  heigh  of  10  meter  m( m/sec );  \  _  /  oun 

radiation  classes. 


Wind  velocities  U.  .  and  t  ,  ,  under  different  wind  directions 

kml  kml 

and  stabilities  can  be  obtained  from  weather  station  routine  statisti¬ 
cal  data.  Wind  velocity  at  stack  height  can  be  derived  using  index 
rate. 


he  formula  employing  ground  routine  meteorological  data  proposed 


by  Nozaki  is  adopted  to  estimate  atmospheric  mixing  depth 


J5] 


Nozak  i  1  s 


formula  is  as  follows: 


where:  T  and  T  are  ground  temperature  and  dewing  point,  respectively, 

U(-)  is  the  average  wind  velocity  at  height  Z,  Zq  is  ground  roughness, 
f  is  the  Clark  coefficient  and  P  is  the  Pasquill  stability  classifica¬ 
tion  .L  from  P=1  for  unstable  class  A  to  P=6  for  stable  class  F). 

Said  formula  includes  contributions  of  thermal  and  dynamic  turbulence. 
The  average  values  of  L  for  neutral  to  unstable  structure  are  obtained 

separately  to  calculated  L  value. 

L  -  m 

III.  Example 

In  order  to  evaluate  possible  pollution  to  the  environment  by- 

sulfur  dioxide_  emitted  from  the  180  meter  stack  at  a  power  plant  which 

is  located  downstream  of  the  Yangtze  River,  the  ground  level  long-term 

average  concentration  was  estimated  using  the  model.  The  terrain 

surrounding  the  plant  is  open  and  flat,  and  is  slightly  inclined  from 

northwest  to  southeast.  The  one-year  routine  data  from  the  local 

weather  station  are  collected  and  statistics  classifications  of 

stability,  wind  velocity,  wind  direction  and  mixing  depth  are  performed 

four  times  a  day  using  the  aforementioned  method  to  obtain  values  of 

parameters  U  ,  <p  and  L  which  are  required  for  the  model  estima- 

tion.  Wind  velocity  is  classified  into  five  levels  based  on  local 

meteorological  features  and  they  are  listed  in  Table  3.  Table  4  gives 

the  statistical  results  of  the  annual  ,  values  under  east  wind. 

kml 

The  average  mixing  depths  under  various  stability  conditions  A,  B,  C, 

D  are  1,400  meters,  1,200  meters,  900  meters  and  650  meters,  respec¬ 
tively. 
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Table  3 . 


Wind  Velocity  Classification  Table 

(  >_ ;  R.iHi  3  & 

(  a  )  i  <  iil  ,'lilli'  •"»/  S  '< 

( 3 )  f  l-J 

Key:  (1)  Wind  velocity  class;  (2)  Wind  velocity  range;  (3)  average 

sind  speed. 


u  1  2  :(  l 

1.0  i.O-l.f  1.5- 4.4  4.5-<.C  C.O 

■ . 0 _ K7 _ 3.5  6.3  7.7  ~  ~ 


Table  4.  Annual  <t>  Values  (east  wind) 

#  mkl 


(i)  e&at 

*•*/ 

) 

A  — D 

c 

0 

* 

C-F 

£ 

1.0. 

0.000‘J 

u 

0.0011 

0.0040 

0.0060 

1.7 

0.0055 

o.oo:u 

0.0076 

0.0172 

0.0337 

5.5 

0.0048 

0.007G 

0.0172 

0.0085 

0.0379 

5.1 

0 

0.0014 

0.0103 

0.0014 

0.0131 

7.7 

0 

0 

0.00C2 

U 

0.0062 

0.0112 

0.0124 

0.0424  { 

0.030s 

0.0969 

Key:  (1)  stability;  (2)  wind  velocity. 


The  plume  rise  and  dispersion  parameter  used  in  estimating  this 
example  are  obtained  from  environmental  wind  tunnel  and  on-site 

■  •  •  r  4 1 

equilibrium  balloon  tests1  J.  The  diffusion  coefficient  adopts  an 
index  form: 


<*  =  C  d 

z  x 


Coefficients  C  and  d  are  constants  related  to  stability  and  are  shown 


in  Table  5 


[4] 


Briggs  formula  series  are  still  used  to  calculate  plume 
ri-se.  Under  neutral  and  unstable  conditions,  j  '  power  is  used  and 

11  i  ii 

the  coefficient  is  1.14;  under  stable  conditions,  j  power  is  used 
and  the  coefficient  is  2.3.  They  are  lower  than  the  1.6  and  2.9  given 
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.  Q  .  [  6  ]  [  7  ] 

by  Briggs 

Table  5.  C,  d  Parameter  Values 

C  D  E  —  F 

0.34 _ -0.-39  -  -  0.16 

0.86  0.80  0.76 

Key:  (1)  Stability;  (2)  C,  d  value. 

The  change  in  wind  velocity  with  respect  to  altitutde  follows  the 
power  rule.  The  on-site  measured  results  using  the  bivane  small  bal¬ 
loon  anemometer  show  that  the  powers  from  unstable  class  A  to  stable 
class  E  -  F  are  0.10,  0.13,  0.18,  0.25  and  0.35  in  that  order. 

The  assessment  area  is  of  a  small  scale  due  to  the  consideration 
of  a  flat  terrain  surrounding  said  plant.  Ignoring  the  effects  of 
terrain,  decay  and  precipitation  scavenging,  the  annual  and  seasonal 
average  concentration  distributions  are  estimated  using  formula  (3) 
for  said  emission  source.  See  Fig.  1.  Said  figures  give  the  con¬ 
centration  distributions  in  an  area  within  a  radius  of  15  kilometer 
with  the  emission  source  as  the  center.  The  low  value  regions  of 
annual  average  concentration  fall  in  the  areas  NE  -  NW  and  S  -  SE  of 
the  source,  and  the  situations  for  spring  and  summer  are  similar  to 
that  of  the  annual.  The  low  value  regions  for  fall  and  winter, 
however,  fall  in  the  area  NNW  -  NE  of  the  source.  The  high  concentra¬ 
tion  value  regions  are  primarily  distributed  in  areas  west  of  the 
source:  both  spring  and  summer  lean  toward  NW,  and  fall  and  winter 

lean  toward  SW.  The  annual  average  high  value  center  is  to  the  NW  of 


C  ,  <l{fc(2) 
C 
d 


A-B 

0.40 

0.91 
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the  source  with  a  weak  high  value  center  located  to  the  east.  The 
distance  from  the  high  value  center  to  the  source  is  shorter  for 
spring  and  summer,  about  2-3  kilometers;  and  it  is  farther  for  fall 
and  winter,  about  7-8  kilometers. 

The  characteristics  of  the  concentration  field  are  closely 
related  to  the  flow  field.  Said  area  has  larger  annual  average  wind 
speeds  in  the  SE  and  NW  directions,  and  they  are  more  evenly  distrib¬ 
uted  in  the  rest  of  the  directions.  Wind  direction  is  predominately 
from  the  east  throughout  the  year  with  southeastern  wind  in  spring 
and  summer  and  northeastern  wind  in  fall  and  winter.  See  Table  6  for 
details.  Figure  2  gives  the  roses  for  annual  wind  direction,  wind 
velocity  and  pollution  coefficient.  From  this,  it  is  obvious  that 
the  concentratT ion  high  value  center  should  be  west  of  the  source. 

The  atmospheric  unstable  centers  for  both  spring  and  summer  are  closer 
to  the  source,  and  they  are  the  opposite  for  fall  and  winter. 


Fig.  1  a . 


Table  6.  Average  Kind  Speed  and 


Kind  Direction 


U )  * 

m>k  Kit  (v: 

IV 

) 

it 

(2) 

(3) 

1 

a  ! 

(4) 

*  t 

£ 

(5) 

if. 

(6) 

(8) 

3.3 

3.1 

2.5 

3.2 

3.0 

NE-Sfc 

(9) 

49 

58 

47 

39 

48 

NE  —  IS 

(i'A 

24 

29 

3G 

31 

30 

E-SE 

&& 

33 

4X 

20 

18 

26 

Key:  (1)  Season;  (2)  Spring;  (3)  Summer;  (4)  Fal^;  <3)  Winter; 
(6)  annual;  (7)  Wind  direction,  Wind  velocity;  (8)  average  wind 
velocity  (m /sec);  (9)  percentage  of  occurrence  frequency. 


Fig.  2.  a.  Annual  average  wind 

velocity  wind  direction 
roses 

Key:  (1)  Wind  speed;  (2)  Wind  direction. 


b.  Annual  pollution 
coefficient  rose 


Although  the  relative  location  of  the  plant  and  residential  areas 
can  be  qualitatively  determined  from  the  wind  direction  and  pollution 
coefficient  roses,  the  distance  from  the  source  and  distribution  range 
of  the  residential  areas  cannot  quantitatively  be  obtained.  The  model 
estimation  considers  the  combined  effects  of  pollution  meteorological 
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conditions  and  dispersion  parameters,  and  therefore,  is  able  to  give 
quantitatively  the  concentration  distributions  of  areas  under 
assessment . 

IV.  Conclusions 

The  method  of  estimating  the  ground  level  long-term  average 
concentration  from  the  power  plant's  high  stack  emission  source  pro¬ 
posed  by  this  paper  has  the  following  advantages: 

1.  The  physical  concepts  of  the  Gaussian  plume  model  are  clean. 

The  calculation  workload  is  not  heavy  so  that  the  capabilities  of  a 
microcomputer  can  be  fully  exploited  Pc— 1500 

2.  It  only  requires  routine  meteorological  data  from  ground 
weather  stations  to  estimate  the  required  pollution  meteorological 
parameters  such  as  wind,  stability,  mixing  depth,  etc.  thereby 
conserving  the  time,  man  power  and  material  required  for  field  observa¬ 
tion. 

3.  It  takes  full  advantages  of  the  existing  research  results  of 
diffusion  coefficients  and  plume  rise  rules.  It  does  not  require  on¬ 
site  field  experiments  of  dispersion  and  plume  rise  in  areas  where  the 
terrain  conditions  are  not  too  complicated. 

The  analytical  results  of  the  example  indicate  that  the  annual 
and  seasonal  long-term  average  concentration  distributions  can  be 
estimated  rather  conveniently  by  the  model.  Based  on  this,  it  will  be 
very  easy  to  conduct  the  plant  site  selection  and  determination  of 
reasonable  distribution  of  various  facilities. 
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C513  ARRADCOM  I 

C535  AVRADCOM/TSARCOM  l 

C539  TRASANA  •  I 

C591  FSTC  4 

C 61 9  MIA  REDSTONE  1 

D008  NISC  1 

E053  HQ  USaF/INET  1 

£404  AEDC/DQF  l 

E406  AFUL  1 

E410  AD/IND  1 

E429  SD/IND  •  1 

P005  DOE/ ISA/DDI  I 

P050  CIA/OCR/ADD/SD  2 

AFIT/LDE  l 

FTD 

CCN  V 

NLA/  PHS  1 

LLNL/Cod*  L-389  l 

NASa/NST-44  i 

NSA/ I213/TDL  : 

aso/ftd/tqia  , 


FTD- ID ( RS)T-0689-86 


